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' Abstract. We study the resonant set of a two-level Schrodinger operator with a linear 

, conical intersection. This model operator can be decomposed into a direct sum of first 
order systems on the real half-line. For these ordinary differential systems we locally 

■ ■ construct exact WKB solutions, which are connected to global solutions, amongst which 

Mh , are resonant states. The main results are a generalized Bohr-Sommerfeld quantization 

' condition and an asymptotic description of the set of resonances as a distorted lattice. 
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1. Introduction 

This paper is devoted to the semiclassical distribution of resonances of the two-dimen- 
sional and two-level Schrodinger operator 



Such two-level operators appear naturally in the study of molecular spectra: if the 
positions of a molecule's nuclei are denoted hy x £ M", and < /i ^ 1 is Planck's constant 
devided by the square-root of the nuclear mass, then a full molecular Hamiltonian reads 
I as H^oi = —h'^Ax + Hei{x). For every nucleonic position x, the electronic Hamiltonian 

■ Hqi{x) is an operator on the electronic degrees of freedom. The full operator //moi acts 

[ on nucleonic and electronic degrees of freedom, that is on wave functions in L^(M^,C), 

^ ' where N = 3 -(number of nuclei-|-electrons) is a notoriously large number. If one considers 

H ■ two eigenvalues of the electronic Hamiltonian Hq\[x), which are well separated from the 

^ I rest of the electronic spectrum a{Hf,\{x)) uniformly for all x, then Born-Oppenheimer 

approximation allows to reduce the study of the full molecular problem to the case of two- 
level systems acting solely on the nucleonic degrees of freedom. The justification of this 
' approximation can be found in |H1 |22j for the time-independent and in (16. for the 

time- dependent case. For considerations from the point of view of theoretical chemistry, 
we also refer to 0. 

Chapter 2 in the monograph gives the standard classification of matrix Schrodinger 
operators with eigenvalue crossings of minimal multiplicity, where our model operator P 
appears as the normal form for a codimension two crossing: the real symmetric potential 
matrix V{x) = V{xi,X2) depends smoothly on the two real parameters xi and X2- It has 
two eigenvalues 

±\J xl + x1 = ±\x\ , 

which coincide on the codimension two manifold {x = 0} C M^. The graph of the mapping 
X I— > ib|x| shows two cones intersecting at the origin, which explains the term conical 
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intersection, which is used in the chemical physics' hterature (see for example |^ llUl05] ). 
Such a degeneracy, geared by two parameters, is generic in the sense, that it cannot be 
removed by symmetry preserving perturbations. The matrix V{x) is in essence Rellich's 
celebrated example of a smooth matrix, which is not smoothly diagonalizable |3Uj . 

If the crossing were not present, one might decouple the system and study the two scalar 
Hamiltonians 

= -h'^Ar, ± \x\ , 

see jS^ISnill]. is a Schrodinger operator with confining potential, and one has pure 
point spectrum only (see Appendix^. The lower level operator P~ , however, has a 
linearly decreasing negative potential, and one sees by a Mourre-type argument, that P~ 
is of purely continuous spectrum. The full operator P inherits the continuous spectrum 
of P~ , while spectrally echoing the discrete spectrum of P~^ with resonances close to the 
real axis. 

In the physical literature, zero energy wave functions close to a conical level crossing 
O El have been studied, and conical intersections have also been addressed as generators 
of resonances jS]. 

A first mathematical proof of existence of resonances for the matrix operator P has 
been given by one of the authors j27j . She crucially used that P is unitarily equivalent to 
the direct sum of ordinary differential operators 

where we use the notation D,. = —id /dr. This decomposition is achieved by a /i-Fourier 
transformation, a change to polar coordinates, a Priifer transformation, and a /i-Fourier 
series ansatz in the angular variable. In particular, r is the length \^,\ ol the dual variables 
^ = (^1,1^2) of {xi,X2) and v is the quantum angular momentum. 

The energy surfaces of the Hamiltonian flow associated with the eigenvalues it 
+ i/2/r2 of the symbol Pu{r, p) of Py are 

{(r,/j) G M+ X M; p2 = {E - r'^ f - jr'^] . 

If < 4E^/27, then the function (E — r^)^ — jr^ has three distinct positive zeros 
< ro < ri < r2, see Appendix 1X1 The Hamiltonian flow consists of two components, a 
periodic flow passing through (ro, 0) and (ri, 0) and a flow coming from infinity and going 
away to infinity passing through (r2,0). Let > be positive. Then rg tends to while 
r\ and r2 tend to \fE as /i ^ 0: the periodic and the unbounded component approach 
each other in the semiclassical limit. 

The resonances of the operator P are defined as the eigenvalues of the complex scaled 
Hamiltonian 

Pq = -hh-^'^^A + e'V(x) , 

which is a non-selfadjoint operator with discrete spectrum independent of the dilation 
parameter 6 G]0, ^[ (see for example 012011^135 ^01 general theory of resonances and 
j27j for the model operator P). In terms of the reduced operators P^, resonances are 
characterized as follows: 

Proposition 1.1 E £ C is a resonance of the operator P = —h^A + V if and only if 
there exist v G /i(N — i) = {|,^,...} and a non-trivial solution w to the equation 

(2) PyW = Ew 
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satisfying 

(3) limu;(r)=0, r^wie-'^r), w' (e-'^r) e L'^{M'^ ,C^), 

for some 9 g]0, 7r/3[. 

The structure of the resonant set can now be analysed by constructing a solution u of 
the undilated system P^u = Eu for E £ C, which vanishes at r = and is incoming at 
r = +00 (see Proposition 14. 

Let Soi{E, h) be the action integral for the periodic flow defined by 



where the square root is defined to be positive when E is positive. It is extended ana- 
lytically into a complex neighborhood of {E > 0}, and if \E\ is bounded from below, the 
asymptotic behavior of Sm{E, h) is given by 

(4) Sm{E,h) = liE^/'^ + ^wh + 0{h'^\\nh\) (/i ^ 0). 

This yields the following Bohr-Sommerfeld type quantization condition for each angular 
momentum P G N — |: 

Theorem 1.2 Let > and z> e N- i = {\, |, • • •} he given. Then there exist e > 0, 
ho> and a function 5{E, h) : {{E, h) e C x ] \E - Eo\ < e, Q < h < Hq] C with 
6{E, h) ^ uniformly in E as h ^ 0, such that equation (0) has a non-trivial solution 
satisfying J^) if and only if (E, h) satisfies the following quantization condition: 

(5) ^ V e-'-/^ E--'/" e2^«i (^'^)/'^ + 1 = 5{E, h) . 



Combining this result with Proposition 11.11 and @, one obtains the following theorem 
about the distribution of resonances. Here, we take A = E'^/'^ as spectral parameter and 
look for resonances in {A G C; a < Re A < 6, ImA < 0, ImA = o(l) as /i ^ 0} for 
arbitrarily fixed positive numbers a, 6 > 0. For each z> G N — |, we define 

^u{h) = |a G C; A = Xkvh - i | {h\n \ - /tin ^x,-}^ , G Z s.t. a < Xkuh < , 

where Xku = ff (8A; - 4z> + 5). 

Theorem 1.3 Let a,b > 0. Then, for any N £ N there exists ho > such that for 
any < h < ho and X G Uf^^N^i^i^) ^^^^^ Is a resonance E of the operator P with 
X _ £;3/2 ^ ^^^j) uniformly for all X G Us<Ar ^u{h) ash^O. 

Remark 1.4 The integer parameters k £ Z should be large of 0{h^^) since a < X^uh. 
Hence, the second term of the imaginary part of A G Ti^{h) is of 0{h) and smaller than the 
first term. Thus, Ti^{h) is an almost horizontal sequence of complex points in the A-plane, 
and UiXAT is a lattice which consists of A'^ horizontal sequences. Theorem II . 31 means 

that for a fixed positive interval [a,b], we can find as many horizontal sequences as we 
want for sufficiently small h, which are close to resonances of the operator P. O 
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Figure 1. Resonances of the operator P = —h'^A + V{x). The parameter 
k lies in {11, 12, ... , 60}, while ly is chosen in {1.5/i, 2.5/i, . . . , 5.5/i}. The 
semiclassical parameter h varies from 10~^ to 1. 

The plots in Figure ^ illustrate the distorted lattice of resonances given by Theorem 
11.31 The larger the angular momentum number u E N — ^, the closer the resonance to 
the real axis. That is, larger angular momentum numbers are associated with longer life 
time of the corresponding resonant states. Studies of the dynamical properties PTIO^I^ 
of the model operator P complement this observation by a quick heuristic argument: 
The symbols p^{x,^) = it |x| of the one-level operators P^ = — /i^A it \x\ induce 
Hamiltonian systems 

\x\ 

with central field. Such systems conserve angular momentum x A ^ = — 2:2^1. 

Hence, only trajectories with small angular momentum come close to the crossing man- 
ifold {x = 0}. Moreover, angular momentum encodes closeness of trajectories to the 
crossing. Straightforward arguments yield that the Hamiltonian system associated with 
p+ has constraint motion including periodic orbits, while the motion corresponding to p~ 
is unbounded. Hence, on a heuristic level, the classical dynamics of the decoupled systems 
reflect the structure of the full operator's set of resonances: a high angular momentum 
number z7 of a resonance mirrors a periodic orbit of the upper level with high angular 
momentum. Such orbits in turn imply existence of localized quasimodes and long-living 
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resonant states. On the other hand, small angular momentum numbers v correspond to 
orbits close to the crossing manifold. Nearby the crossing, non-adiabatic transitions to 
the unbounded motion of the minus-system are possible. Hence, in this regime shorter 
life-times and resonances far away from the real axis have to be expected. This heuris- 
tic point of view is in wonderful agreement with the derivation of the Bohr-Sommerfeld 
conditions ©, and the prefactor 



can be seen as some type of Landau-Zener probabilty for not performing a non-adiabatic 
transition near the crossing. Hence, the conditions © is an interpretable extension of 
quantization conditions for scalar Schrodinger operators to the case of a matrix-valued 
operator with crossing eigenvalues. For a study relating the resonant set of the matrix 
operator P with the spectrum of the scalar operator we refer to Appendix IbI 

The method of proving Theorem 11.21 and 11.31 is similar to jl3j , while the problem re- 
sembles in some part the double well problem We construct exact solutions of the 
equations PyU = Eu, v € /i(N — |), which are exact WKB solutions in local domains 
of the complex plane. Connecting these solutions from domain to domain, we get global 
solutions. More precisely, we construct an exact solution which vanishes at the origin, 
and represent it, after several connection procedures, as a linear combination of Jost solu- 
tions, which are defined at infinity. The quantization condition of resonances will be given 
as the condition that the connection coefficient of the outgoing Jost solution vanishes 
(Proposition 14. 1|) . 

The proof of our main results. Theorem 11.21 and Theorem 1 1.3( involves several different 
tools and proceeds as follows: In ^ we reduce the study of the full operator P to that 
of the ordinary differential operators Py^ v G /i(N — ^), and prove Proposition 11.11 In 
331 we extend the theory of exact WKB analysis for one-dimensional scalar Schrodinger 
operators used in jl5) to a class of Schrodinger systems. This exact WKB theory works 
in the generic situation, where the operator is without singularity, and turning points are 
"sufficiently well-separated" from each other with respect to h. Our systems, however, 
have the origin as a regular singular point. Moreover, the first turning point ro tends to 
this singularity, while the other two turning points r\ and r2 approach each other at \fE 
as /i ^ 0. ^outlines the strategy to obtain a global solution under these circumstances. 
In ^ we define Jost solutions and represent them as exact WKB solutions. In ^ we 
construct an exact WKB solution vanishing at the origin, facing the same kind of difficulty 
as in Jl]. The bad error estimates in Theorem II. 21 and Theorem 11.31 come only from here. 
33 is devoted to the connection formula at \fE. Here, we reduce the operator Py to a well- 
studied microlocal normal form. In 3E1 we compute the quantization condition (Theorem 
II. 2|) from the connection formulae obtained in the preceding sections. This condition 
is given in the form of Bohr-Sommerfeld using the action Soi{E,h) between ro and ri. 
Analysing the asymptotic behavior of the action Soi{E, h), we finally get the semiclassical 
distribution of resonances (Theorem II. 3|) . 

Acknowledgments. Part of this work has been carried out while the authors have 
been members of the program on semiclassical analysis at the MSRI, Berkeley. We thank 
Maciej Zworski for valuable discussions. 
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2. Reduction to the first order ordinary differential system 



The aim of this section is twofold: first, we define the resonances of 

(6) P = -h^A + V{x), V{x)=(^ll 

as the eigenvalues of its complex scaled counterpart 

Pg = -e-^"^h^A + e''^V{x), ^G]0,f[. 

Then, we reduce the study of the two-dimensional operators P and Pg to countably many 
one-dimensional operators {Pu}i/eh{i/2+z)j 

— hDr u/r 
v/r r"^ + hDr 



Pu 

and {Pe,u}ueh{i/2+z) 



where stands for —id /dr. Most of the material presented here can be found in |27j . 
We give a self-contained account of it for the convenience of the reader. 

Lemma 2.1 Let 9 g]0, The operator Pq with domain 

V={ue H'^{R'^, C^) ; \x\u G L'^{R^, C^)} 

is closed and has purely discrete spectrum, a^Pg) = o"disc(-P6»)- The spectrum is indepen- 
dent of the scaling parameter 6 g]0, ^[ in the sense that 

(Tdisc(^'e) = (TdM') {0<9<e' <^). 



Proof: We split the proof in three steps. 
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Closedness. Let a G C with Ima ^ and Tq, := — /i^A + aV{x). Then, the following 
quadratic estimate holds: there exist positive constants Ca,ba> such that 

(7) \\To.u\\h + Ca\\u\\l2 > ba (II - h^^u\\l2 + II \x\u\\l2) 

for all u E C^{M?,C'^)- Indeed, denoting D = -iV^, 

TaTa = {hDf + Im(a) V{hD) + Re(a) (-/i2Ay(x) - V{x)h'^A) + |ap|xp , 

where one uses i[— /i^A, = V{hD). Since = |xp, one gets 

< (-/iVr^^ ± |a| = |ar^(/iZ))^ ± (-/i^A)y(x) ± y(x)(-/i2A) + |a| |xp 

and 

Re(a) (-/i2Ay(x) - V{x)h'^^) > -\ Re(a)| \a\~^ {{hOf + |Qp|x|2) . 
Since V{-) behaves like ±| • |, one obtains 

To^T^ > {l-\Re{a)\\a\-'^){{hD)'^ + \af\xf) +lm{a)V{hD) 
> {{hDf + |x|2) - c„ 

and the desired estimate ((T)). Literally the same arguments as in the proof of Theorem 
IL3 in 20 conclude the proof, that Ta is a closed operator on V. 

Discrete spectrum. Let 6 g]0, -^l. Lemma 3.1(ii) in |27j proves, that the mapping 
Pe - ze -.V ^ L2(M2^C2) is bijective for ah G C with Im(zee-*^) > 0. By Rellich's 
compactness theorem, the embedding V — > L^(]R^,C^) is compact. Hence, {Pq — zq)'^ is 
compact. Let S C with Im(zee~*^) > 0. For z S C one writes 

{Pg - z){Pe - Zg)-^ =I+{ze- z){Pg - Zg)'^ =: I + Kg{z) , 

where {Kq{z)} z&c is an analytic family of compact operators. For z sufficiently close to 
zg, one has ||i^e(-z)|| < 1- Then, by the analytic Fredholm theorem, z ^ {I + Kg{z))~^ 
and z I— >• i^Pg — z)^^ are meromorphic in C, and the residues at the poles are finite rank 
operators. Hence, (y{Pg) = Cdisc(-Pe)- 

Independence of 6. For G S = {0 G C; lm.4> G]0, |[} one defines 

:= -e-^'^h'^A + e^V{x) 

as a closed operator with domain V. Since Hfj, = U^^^Pg U^-iy.^ with 9 = Im cp and 
URe<i>u{x) = e^'^'l'u{e^'^'^x) the unitary scaling, one has cf{H^) = a{Pg). Since {H^^^^^s 
is a holomorphic family of type (A), eigenvalues of depend analytically on cj), see 
VH.§3 in By unitarity, (j) ^ is constant if lm.(j) is constant. Hence, the eigenvalues 
of H(j, and Pg are independent of (j) and 0, respectively. □ 

We note, that the constants 6q,Cq > used in the preceding proof approach zero as 
Ima tends to zero. Hence, the above quadratic estimate does not yield closedness of the 
undilated operator P = —h'^A + V{x) on the domain V. However, a straight forward 
adaption of the Faris-Lavine Theorem proves essential self-adjointness of P on D. 

Definition 2.2 The eigenvalues of the dilated operator Pg = — e~^'^/i^A + e'^^V{x), 
6 g]0, f[, with domain P = {« G ^^(M^C^); |x|n G L'^{R^,<C?)] are called the resonances 
of the operator P = -h?A + V{x). 
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The choice of dilation angle 9 g]0, ^[ is crucial: under x i— > —x the dilated operator 
= — e^^*'^/^/i^A + e'^^^^V{x) is unitarily equivalent to the original operator P, which 
has no discrete spectrum, see Proposition 15.21 

Lemma 2.3 Let us fix 9 €]0, Then E G C is a resonance of P if and only if there 
exists V E + Z) such that E G C is an eigenvalue of the operator Pq ^ with domain 

if V ^ I , and 

V^H ={we L2(M+,c2); {-DrWi ± {2r)-^W2), {DrW2 ± (2r)-^u;i), r'^w G L^{R+,C^)] . 



Proof: Again we proceed in several steps. 

Fourier transformation and polar coordinates. Let n(^) = {J-hu){^) be the semiclassical 
Fourier transform of u, 

n(e) = I e-'''<l^u{x)dx. 

Then the dilated equation 

(8) Peue = Eue 

becomes 

Switching to polar coordinates (■^1,^2) = r(cos sin 0), r G M+, (/> G T = R/2tt'L one gets 
|e-2*^r2 - hA{(t))e'^Dr - ^^{(p)e'^D^ - e| = 

with 

cos (b sin 6 



A{<t)) -- 

' \ sm.<p — cos( 

Since the volume element dx changes from d^ to rdrdcj), the new Hilbert space is 
{u G L2(M+ X T,C2); G L2(M+ x T,C2); ^(r, (/> + 27r) = u{r,(p)} . 

Half-angle rotation and conjugation by \pr . Now put 

WQ = y/rR{-^)u0, 

where 

Y sm cos (, 
is the rotation by the angle (p. Using the relations 

A{<P) = R{<P)i^ \ -l) = (o i?'(<A) = i?(<A+f), 

we obtain the equivalent equation to © as 

(9) I _ ^..^ ^1 /) z)^, _ ,^oh / M ^ ^ = 0. 



^ -1 y ' r V 1 . 

The new Hilbert space is [w G L2(M+ X T,C2); w{r,(t) + 2tt) = -w{r,4))]. 
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Fourier series expansion. We expand w$ in a Fourier series, taking into account that 
wg is 27r anti-periodic with respect to (j): 



Then we see that the original eigenvalue problem Q is reduced to that of the differential 
expression 

The resulting Hilbert space is L^(M+,C^). 

Closedness of Pq ^ on V^. As before in the proof of Lemma l2.1l one shows a quadratic 
estimate, this time for the operator 



2 



with a € C, Im(a) ^ 0: there exist positive constants ba,iy,Ca,u > such that for all 
u;GC~(M+,C2) 

b^^^ {WhDrwWl^ + \\r-'^w\\l2 + Wr^wWl^) if / ±| , 



\Ta,uW\\h +Ca,u\\M\h ^ ^ 



-hDr - \ _ _j^2 



r 



3h I 
2 ' 


5h 
2 ' 


u 




hDr 










± 


hDr 





Ar + ^ + ^h^(^^ > -/i^A, + ^ , 



and 



Since v'^ — h\u\ > for u ^ ib|, one obtains for suitable ^, Ca,u > 

r^,.T„,, > (l-|Re(a)||ari) (-/i2A, + ^^^ + |a|V) -|Im(a)|2r/i 
> ba,,y {-h^Ar + r'^ + r^) - c„,^ . 
The estimate for u = ziz^ follows analogously. □ 

Proof of Prop osi tion U.ll Let be a resonance of P = -h^A + V. By Lemma ITTH there 
exist u £ h{Z + 1/2) and wg £ Vy such that Pe,yWe = Ewg. The origin r = is a regular 
singular point for the operator Pg^i, with indicial roots ±9 = zizv/h. Hence, by the theory 
of Fuchs, w$ is a linear combination 

we{r) = Cow'>{r) + CooW^{r) 

of the two solutions w^{r) and w^{r) to the equation Pe^yW = Ew such that 

w\r)^r\''\( .y w^{r)^r~\''\( ]...\ (r ^ 0), 

^ ' \ — sgn(z^) ^ y ' ^ ' y sgn{u)i J ^ " 
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see also §4. Then, the condition wg G imphes Coo = 0, and we behaves hke r''^' near 
the origin. In §5, the construction of complex WKB solutions with base points at infinity 
shows 

where w^{r) is exponentially growing and w^i^r) exponentially decaying as r ^ oo. 
Hence, C+ = 0, and we{r) decays exponentially as r — > oo. Setting w{r) := wg{e'^^r), one 
obtains a solution to P^w = Ew with the claimed properties. 

Let E £ C and w he a solution of PuW = Ew for some u G h{'L + 1/2) such that 
\mir^o+w{r) = and r^z/;(e~*^r), tt;'(e~*^r) G L^(M"'",C^). The preceding arguments 
yield, that we{r) := w{e~^^r) is a solution to the equation Pe^yWe = Ewq with we G Vy. 
By Lemma 12.31 ii^ is a resonance of P = —h?IS. + V . □ 

Remark 2.4 Let E G C and v G h{'L + ^). Then the following equivalence holds: 
u = (ui,U2) is a solution of P^u = Eu, if and only if -u = (— ui,U2) is a solution of 
P-uU = Eu. The same holds for the dilated operators P^fi and P~u,e, ^ £]0, Hence, 
we will restrict our studies to the case i/G/i(N — ^). O 



3. Exact WKB method for 2x2 systems 

We now wish to find a representation formula for the solutions of Py, from which it is 
possible to deduce the asymptotic expansion in h. The method is known as exact WKB 
method. We derive it in a somewhat more general context, and then apply it to our 
specific equation. 

We study 2x2 systems of first order differential equations in a complex domain 
which are of the form 



or equivalently 

I dx \ -uj{x) -p2{x) J 

The functions pi,P2, and uj are holomorphic in D. The following considerations will lead 
to the construction of exact WKB solutions for this type of systems. 

3.1. Formal construction. A usual change of variables (using oscillatory part) and some 
other basic transformations reduce the operator to a more computable one. Let us put 

g+{x) = |(pi(x) +P2{x)) +uj{x), 5_(x) = -i(pi(x) +P2{x)) +uj{x). 

After conjugation by 

M{x) = exp r{pi{t)-p2{t))dt] f \ ]]=: mix) ^ ^ ^ 



the system (|10j) is transformed into the trace-free system 
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with u{x) = M{x)v{x). Introducing a new complex coordinate 

(12) z{x) = z{x;xo) = V9+{t)9-{t) dt , xq e D . 



XO 



we look for solutions of the form e hw±{z). 

Definition 3.1 Let pi, p2, and uj be holomorphic functions in D. The zeros of the 
function g^{x)g-{x) = —\{pi +^2)^ are called the turning points of the system (fTII]). 



We note that due to the possible presence of such turning points the square root in 
the definition of z{x) might be defined only locally. By formal calculations, the amplitude 
vector 'w±{z) has to satisfy 

h d ^ ( ±i H{z)-^ \ ~ ( ^ 

where the function H is given by 

For a decomposition with respect to image and kernel of the preceding system's matrix, 
we conjugate by 

p _ 0-1 ( H{z) ±iH{z)-' \ 1 _ / H{z)-^ H{z)-^ 

^^^^^ ~ ^ \^ H{z) ^iH{z)-^ y ' ^± - y ^iH{z) ±iH{z) 

and obtain a system for w±{z) = P±{z)w±{z), 

a , . / u Hiiz) 



where H'{z) is shorthand for dH{z) / dz. The series ansatz 
(13) w^{z) = Y, 



W2n,±{z) 
W2n+l,±{z} 

n>0 

with Wo, it = 1 and for n > 1, the recurrence equations 

(14) U^7lj"-+^'^(^) = ^"-'^(^)' 

(15) J^W2n+2,±{^) = W2n+l,±{^) 

give us a formal solution up to some additive constants, which are fixed by setting 

Wn,±(z) =0, n > 1, 

for a base point z = z{x) where 5; E D is not a turning point. We note that the preceding 
equations for Wn,± are the same as the ones obtained by an exact WKB construction for 
scalar Schrodinger equations. See for example the work of C. Gerard and A. Grigis 
or T. Ramond 29 . 

Let be a simply connected subset of D which does not contain any turning point. 
Then the function z = z{x) is conformal from O onto z(0). Assume that z E z(0). If 
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r±(z, z) denotes a path of finite length in z{0,) connecting z and z G z{Q), we can formally 
rewrite the above differential equations for n > as 

f H'(C) 

W2n+l,±(.z) = / exp{±l{C - z)) ——W2n,±{C)dC, 

Jr±{z,z) ^lU 

W2n+2,±{z) = / T777T^2n+l,±(C)rfC 



or after iterated integration as 

^^2n+l,±(^) = 



r±(z,z) 7r±(?,C2n+i) 



/ exp(±f(C2-C3 + --- + C2n+l-^)) X 

-'r±(2,Ci) 



H{Cl) i?(C2n+l) 



d(i . . . dC2n+l, 



W2n+2,±i^) = / ••• / exp(±|(C2-C3 + ----C2n+2)) X 

ir±(i,^) Jr±(?,C2n+2) ^r±(?,Ci) 

3.2. Convergence, /i-dependence, and Wronskians. We now give the preceding for- 
mal construction some mathematical meaning in turning point-free compact sets Q C D. 

Lemma 3.2 For any fixed h > 0, the formal series i)Jc^|) converges uniformly in any 
compact subset of Q,, and 

(16) W'^'^{x,h) =Y,W2n,±{<x)), Wl^''{x,h) = Y,W2n+l,±iz{x)) 

n>0 n>0 

are holomorphic functions in D. 



Proof: In 17, all the functions defined above are well-defined analytic functions. For 
compact subsets K <Z ^ and J, z G z{K) there exist positive constants C^{K) > 
depending on the semiclassical parameter h and the compact K such that 



(17) 



sup 

C6r±(?,2) 



exp (±10 



H'iC) 
H{C) 



< Cl{K). 



If we denote the maximal length of the paths r-|-(z,-) C in the preceding iterated 
integrations by < L < cxd, then 



sup \Wn,±{z)\ < 
z&z{K) 



n>0, 



where the bound ^ comes from the volume of a simplex with length L. □ 
Thus, we have uniform convergence of the series (|13|) for w±{z) and exact solutions 



u±{x) 



h m{x)T±{z{x)) 



(x) 
w'^'^ix) 



y^even 
,,odd 
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of the original problem ()1U() on turning point free sets O, where 

1 1 \ / H{z)-^ H{z)-^ 



^^^^^ ~ » -1 1 A ^iH{z) ±iH{z) 



(18) 

H{z)-^^iH{z) H{z)-^±iH{z) , 
-H{z)-^^iH{z) -H{z)-^ ±iH{z) ) ' zi,z[iL). 

We write these solutions u±{x) as 

u±{x;xo,x) 

indicating the particular choice of the phase base point xq in (|12j) . which defines the phase 
function z{x) = z{x;xo), and the choice of the amplitude base point z = z{x), which is 
the initial point of the path T±(z, •). 

For a fixed x G let il,± be the set of x G such that there exists a path T±{z(x), z{x)) 
along which x i— > ibRe2:(x) increases strictly. Then, 

Proposition 3.3 The identities for w^^'^{x,h) and w'^'^{x,h) give asymptotic 
expansions in Q± . More precisely, we have for any a G N and S N 

N 

n=0 
N 

n=0 

uniformly in compact subsets ofQ±. In particular, 

<^^'^(x, h) = l + 0{h), wf'^ix, h) = 0{h). 



The proof is just the same as that of Proposition 1.2 of jl5j . The key point is the follow- 
ing: Since the iterated integrations defining Wn,±{z) contain terms of the form exp(ibC//i), 
one has to make sure that C, i— > =bRe(C) is a strictly increasing function along the path 
'^±(z,z). In other words, the paths V ±{z{lic) , z{x)) have to intersect the Stokes lines, that 
is the level curves of x i— > Re{z{x)), transversally in a suitable direction. 

One defines the Wronskian of two C^-valued functions u, v as W{u, v) = uiV2 — U2V1. If 
z = au + Pv with a, /3 S C, then 

_ W(z,t;) W(z,n) 
" W{u,v)' ' W{u,v)' 

Elementary computations give the following exact Wronskian formulas for exact WKB 
solutions with different phase and amplitude base points, using w'^^^ and w'^^. 

Lemma 3.4 Let x, xo,yo,x,y £ $7. Then, 

(19) W{u±{x;xo,x),u±{x;yo,y)) = ±2im{xf exp (±i (z(x; xq) + ^(x; yo))) 

X (<-(x;xo,i)«;5:''''(^;yo,y)-«^±^(x;xo,£)«;r'^(x;yo,y)) , 

(20) W{u±{x;xo,x),u^{x;yo,y)) = ±2im{xf exp {±j^{z{x]Xo) - z{x;yo))) 

X (wT'^ix; xo, SX^'^(x; yo,y) - w°^^\x; xq, iX^(x; yo, y)) . 
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In particular, if pi = P2, then m = 1, all the Wronskians are independent of x, and we 
have for solutions with the same phase base point 

(21) W(u±(- ;xo,x),u±(- ;xo,y)) = exp {±^z{y]Xo)) Wj^'^{y;xo,x) , 

(22) Wiu±{-;xo,x),u^{-;xo,y)) = ±2i w^^y ; xo,x) . 



4. Formulation of the resonance condition 
We now go back to our system 

(23) {P,-E)u{r) = ( ''^-^-^^ .2_i^AAVw = 0' 

V r ' ^ ^ i dr J 

for £^ E C and u S h{H— ^), and study in particular their behaviour as r ^ and r — > oo. 
In what follows, we use x instead of r and rewrite the equation H23|) in the form 

f^A\ , A A ( x^ — E vlx 

(24) hD^u = Au, A=^ 

The origin is a regular singular point of the equation with indicial roots ibz>, where 
z> := > is positive, see Remark 12.41 Indeed, H24() can be rewritten as 

x^u={Aq + 0{x))u (x^O), Ao=(^_°.^ 

and ibP are the eigenvalues of ^-nd *(1,=fO the corresponding eigenvectors. Let 
uq{x) be a solution corresponding to the index v, 

uq{x) ^x"" 

as X ^ 0. We will see in the next section, that there exist Jost solutions /^(x), which are 
characterised by their asymptotic behaviour at infinity, 

(25) /+(x) ~ e-^i{-'-^E.)lzh 1 ^ ^ _ g-.(x3-3i^x)/3h 

as X ^ +CXD. If 9 g]0, 7r/3[, then /+(xe^*^) is exponentially growing and /^(xe~*^) 
exponentially decaying as x ^ +cxd. Since and /" are linearly independent, uq can be 
expressed as a linear combination of these solutions, 

(26) uo{x) = c+{E,h)f+{x)+c-{E,h)r{x). 

^From Definition 12.21 and Proposition 11.11 we obtain the following characterization of 
resonances: 

Proposition 4.1 The energy E £ C is a resonance of P if and only if there exists 
ueh{N-^) with c-^{E, h) = 0. 

To calculate the coefficients c^{E, h), which connect the solution uq defined at the origin 
with the Jost solutions defined at infinity, we need some intermediate solutions, which 
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we will construct as exact WKB solutions. Let us recall the WKB construction of 
this case. Note that tr^ = and so m = 1. Exact WKB solutions are of the form 

(27) n±(x; xo, x) = e±^(^')/'^r±(z(x)) ( ""jZ^^}^ ) , 

where the phase function z{x) is defined by 

^ 

z{x) = z{x; xo) = y/g+{t)g_{t)dt, g±{x) = - ^ E ± x'^ 

for a phase base point xq, T±{z(x)) is a 2 x 2 matrix defined by (jlHj) with 

'5_(x)\^^^ f V + Ex — x^\^^^ 



H(z(x)) — , , , , — , ^ „ 
^ ^ \g+{x) J \i^-Ex + x^^ 

and the symbols w^'"^{x,h) = En>o '^2n,±(2:(a^)) and w^'^{x,h) = Yln>o'^'iri+i,±{z{x)) 
are constructed by recursive integrations with an amplitude base point x. 

There are at most six turning points, the zeros of x i— > g^{x)g- (x), in the whole complex 
plane C^;. They are point-symmetric with respect to the origin and denoted by {i^jj^^Q, 
see Appendix For E > fixed and v > sufficiently small, they are real and satisfy 

< ro < n < Ve < r2. 

As /i — > 0, that is ^ 0, tq tends to and ri and r2 tend to ^/E, where ^/E denotes the 
square root of G C in the right half-plane. Notice that rg and ri are zeros of g+{x) and 
r2 is a zero of g-{x). 

The Stokes curves are the level curves of the function x i-^ Re2;(x). The Stokes curves 
emanating from the turning points are drawn in Figure 2. 




Figure 2. The Stokes curves, that is the level curves of x Re2;(x), 
which emanate from the turning points ro, ri, r2- The arrows indicate the 
directions along which Re z{x) increases. The dashed lines show the branch 
cuts, xf and xf are the amplitude base points for the constructed exact 
WKB solutions. 

We put branch cuts as drawn in Figure 2 for the multi-valued functions 

V9+{x)9-[x) = , H{z{x)) = [ , 

X — Ex + X- 



16 



S. FUJIIE, C. LASSER, AND L. NEDELEC 



and suppose that 



Q \ 1/4 

u + Ex - x'^\ 



V.'-XHE-X^)%=0 = (^ ,_j,^ + ^3 J U=o = 1. 

Let -E > be positive, and let = P/i > be sufficiently small. Then, for example, 
\J gj^{x)g- [x) S iM"*" for tq < x < r\ and ri < x, and H{z{x)) G e~*'^/^ M"^ for ro < x < ri, 
while i^(2(x)) G e*''/^M+ for r2 < x. 

Choosing amplitude base points xf, x^ as in Figure 2 (x^ are supposed to be purely 
imaginary for a technical reason in the proof of Lemma ^ , we define the following six 
exact WKB solutions 

Uq{x) = u±{x;ro,xf), 
uf{x) = u±{x;ri,xf), 
n^(x) = u±{x;r2,x^) 

with the turning points tq, ri, r2 as phase base points. The three pairs {u'^ (x) , Uq (x)) , 
(x)), («^(x), u~ (x)) are all linearly independent and they are connected with 
uq and {f'^,f~) by transfer matrices ^{cq {E,h), Cq {E,h)), Ti{E,h), T2{E,h), T2,{E,h): 

(28) .xo(x) = «(x),no(x))( IllJ^j 

(29) (n+(x),no (x)) = {ul {x),ul {x))Ti{E,h) 

(30) (n^(x),n^~(x)) = (n+(x), (x))r2(i?, /i) 

(31) (n+(x),n;(x)) = (/+(x), /- (x))r3(i?, /i) 
Then, the coefficients c'^{E,h), c^{E,h) in are given by 

(32) (;;!)=T3r,r.(| 

Now we can describe the strategy for the remainder of this article: The matrix T3 for 
the transfer at infinity is computed in Section |S1 the connection coefficients at the 
origin in Section |HJ and the transfer matrix T2 near x = in Section [7| The matrix Ti 
is the easiest one, and can be determined right away using Lemma 13.41 Because Uq and 
uf differ only in the base point of the phase, one uses (fTOj) and (|2fl]) for x = x^ to obtain 



(33) 
where 



Soi{E,h) = / y/g+{x)g^{x)dx = / ^ ^ Ux 



is the action integral between the turning points ro and r\. 
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5. JOST SOLUTIONS 

The Jost solutions of Pi,u = Eu are characterised by their behaviour H25|) at infinity. 
They can be expressed as exact WKB solutions with the base points of both phase and 
amplitude placed at infinity. This fact allows us to calculate Ts. 

First we define the phase function with base point at infinity, 

z(x,oo) = - t-^{E - t2)2/t _ i(f _ E)) dt + - ZEx). 

Taking the branch of the square root into account, we see that the integral converges 
absolutely, hence 

i o 

z{x,co) = -{x — 2>Ex) + o{\) {x^+oo). 

Next we define amplitudes based at infinity. The Stokes curves are asymptotically like 
horizontal lines {Imx = const.}, and Imx i-^ Yiez{x) is a decreasing function (see Figure 2). 
As in Section 3 of we choose infinite paths 'y±{x) starting from infinity and ending 
at X, which are asymptotically like lines of the form {Imx = ^5^ex} for some 5 > 0, 
such that X ^ ^^ez[x) are strictly increasing functions along 7±(x). Denoting the path 
z{^±{x)) by T±[+(X), z{x)) and setting t«o,± = !> we inductively define Wn,±{z) by 



W2n+l,±{^) = / exp{±j^{C - z)) ——W2n,±{0dC, 

Jr±{+oo,z) -HlU 

W2n+2,±{^) = / W2n+l,±iC) dC , n>0. 

Jr±(+oo,z) ^{(,) 



Noticing that 



'r±(+oo,z) 

f 

'r±{+oo,z) 
H'{x) V E-'ix^ 



0{x- 



H{x) 2 {v-Ex + x3)2 

as X — > +00, one constructs well-defined exact WKB solutions u^{x) corresponding to 
these base points, proceeding as in Section 3 of [5^]. Up to a constant prefactor, u^{x) 
are the previously defined Jost solutions: 

Lemma 5.1 Let he the exact WKB solutions with phase and amphtude base point 
at infinity, the Jost solutions defined in i)25|) . Then, 

f^{x) = ±^e^'/^ut{x), x>0. 



Proof: We just check the asymptotic behavior of at infinity. Since H{z{x)) — > e^"^/^ 
as X ^ +0O, we get by an elementary calculation 

1 1 W e-'^*/^ \ / 1 

-1 ill ^ie^i/^ ±ie^^l^ Mo 



(x) ~ g±i(a;^-3£;x)/3h 

that is, 



ni(x) ~ 2e-'^^/^e+^(-^-3^^)/3'^ ( J ) , n- (x) ~ -le 



n/i^-i{x^-1.Ex)/'ih ( 
1 



as X ^ +00. □ 
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As an immediate consequence of the behaviour of the solutions u^{x) for x — > +00, we 
obtain that the discrete spectrum of the full operator P is empty. 

Proposition 5.2 Let P he the full operator as defined in (0). Then, (Tdisc(-P) = O- 



Proof: The proof of Lemma 12.31 shows, that £J G R is an eigenvalue of the self-adjoint 
operator P if and only if there exists v G h{N — i) such that £^ E M is an eigenvalue of the 
ordinary differential operator P,^. At infinity, every distributional solution of (j2SI) has to 
be a linear combination of However, for G M neither /"*" nor /~ are in L^(]R^,C^). 
Hence, the discrete spectrum is empty. □ 

The main result of this section is the following proposition. 

Proposition 5.3 There exists a positive 5 > independent from E G C and h > 0, 
such that the transfer matrix T3 defined in l^31\) satisfies 



T3{E,h) = 2e-^'/^ 



e^2-(^''^)/'*(l + 0(/i)) 0(6-^/^^) 

0{e-^/h) e-^2^(-^''^)/'^(l + 0(/i)) 

as h ^ 0, where S2oo{E, h) is the action between r2 and +00, 



S2ooiE, h) = j^^ (^^u^-x^{E-x^Y/x - i[x^ - ^)) dx + i(r3 - SSrz). 



Proof: We use the Wronskian formulas of Lemma 13.41 and the asymptotic expansions of 
Proposition 13.31 to calculate Ta. Setting 

(n+,^i-) =: (u+,n^)r3, 

the previous Lemma l5 . 1 1 gives 



2g-7ri/4 

-2e-'^'/'^ 



The difference of the phases z{x;r2) and +00) is the action S2ooiE, h), and we have 

Wint,n^) = ±2i e^'^'-'^ (1 + 0{h)) , >V(n± , n^) = ±2i (1 + 0(h)) . 
Since there exists 5 > with 

one obtains 



□ 
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6. ASYMPTOTICS OF THE SUBDOMINANT SOLUTION NEAR THE ORIGIN 

We recall that the origin is a regular singular point of with Fuchs indices itP, and 
that no is a solution corresponding to which is unique up to a constant multiple. The 
purpose of this section is to calculate the connection coefficients c^{E,h) in (EHJ, i-e. to 
connect uq with the exact WKB solutions Uq{x) and ttg (x). For this, we first need to 
express uq as an exact WKB solution. 

Let us look at the asymptotic behaviour of the phase z{x) = z(x;ro) as x ^ for a 
fixed h > 0. Since \J v'^ — x'^(E — x^)^ = + x<f)(x) with 4>(x) holomorphic near x = 0, 



e 



^W/'^ = exp(/ {u/t + ^{t))dt/h] =CE,hx'^exp(^J ^{t)dt/h 



where CE,h ■= r^"^ eyip{- (j){t)dt/h) > is a positive constant. Moreover H(x) is a 
holomorphic function at x = with H{0) = 1. Thus, 

(34) e<^y^T+{z{x)) ( J ) ~ CE,hx' ( J ) , x ^ 0, 

while h > fixed. This suggests that uq is collinear to an exact WKB solution of the 
type +. However, it is important to notice that the amplitude base point x of the WKB 
solution should be placed at the origin in order to have 

n;r"(x) l\ ^ ^ 

w%^^ix) J V / 

since x i— > Re 2:(x) is decreasing as x tends to in radial directions (see Figure 2). Moreover, 
the origin is a singular point for the equation and we need to check that the exact WKB 
solution n_|_(x;ro,0) is well defined, i.e. the recurrence equations ()14|). H15|) with initial 
value t^o ^ = 1, ^(0) = 0, n > 1, define a sequence of holomorphic functions {w^ _^_}n>o, 
and the series X^n>o ^2n,+(^) ^^'^ Sn>o '"^2n-i-i,-i-(^) converge in a neighborhood of the 
origin. We rewrite (|14() . (jlSj) with respect to x: 

(35) + ^V<7+ <+i,+ (^) = 

(36) ^^°"+2,+ (^) = f^^2Vi,+ (^), 

where H!j. stands for the derivative of H with respect to x. 
These equations are of the form 

dw b(x) „, , 

^ + V^ = /(^)' ^(o) = o 

with b{x) and /(x) given holomorphic functions at the origin. In our case, 6 = for H36|) 
or 6(x) = 2xy^ g+ {x)g- (x) / /i, 6(0) = 2i> for (j3SJ). This Cauchy problem has a unique 
holomorphic solution if Re 6(0) > —1, and the solution is given by 

(37) w{x) = xj^ t^(°)exp(^-x^ b{xs)ds^ f{xt) dt, 

where 6(x) = (6(x) — 6(0)) /x. 

Hence for D > 0, are uniquely determined and given by the recursive integrals 

^^2n+2,+ = Ioi'W2n+l,+), W^n+1,+ = h{w2n,+), 
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where 

hif) = exp(^-- V5+(05-(i)^iJ -jf^fiOd^ 
It is not difficult to see that for any fixed positive /i > 0, the series 

oo oo 
?i=0 n=0 

are absolutely convergent in a sufficiently small neighborhood of the origin. Hence, the 
function 



w 





odd 



:o'(^) 



defines a solution to (|24|) . 

Next, we study the asymptotic behaviour of the connection coefficients c^{E,h) as 
/i — > 0, using the exact WKB solution uq. For that purpose we will need some bounds on 
w'^^(^{x, h) and w°^q{x, h). Since we will finally use those in the Wronskian formulas, it is 
enough to deal with the case when x is purely imaginary, x = iR with i? > 0. We start 
by the following elementary estimate. 

Lemma 6.1 For any «; > 0, m > and r > 0, one has 



e 



'('■"^^T ^^dr < 







(1 + r)"^+i ~ m \1 + T 



Proof: One integrates by parts, 

m-l m \ r d f p'^(''"^) \ 

e«('^-^) 1 = - — / r*"— — - \ dr 

(l+r)"^+i m(l + T)™+i m Jo drl(l+r)™+M 

1 / r 

< 







m \1 + 

using that the subtracted integral is positive. □ 
Lemma 6.2 For any E > 0, t > 0, and n G N, one has 

(39) l<+^^'^"/^)l ^ ^ (tTvJ ' 

wiiere i^(r) = 1 + Su'^t'^/E^. 

Proof: For x = iR, R > 0, we have by the changes of variables t = is, ^ = 

2? Jo + p^(F + 
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r,n, 1 f"' / 2 f' ^^^' + a'jE + sY , \ •'(E + V-) , 



1/2 +p2(^ + p2)2 - 1.2 + ^2^2 ' 

we have for R = / E, 

(40) \W)\.=^ur/E\ < ^ r 1/ (^^r/ii;) | dr, 

(41) |/i(/)U=Wi^| <^ j\^'^^-^^r^\f{ivT/E)\dr. 

We now proceed by induction over n G N. Since + = 1, inequality ((HHl) is trivially 
satisfied for n = 0. Next, we assume that (|39|) holds for n = 2k. Then, from (|41|) and 
Lemma l6.H one has 

2k\ Jo (i + r)2fc+2" - {2k + iy. \l + T^ 

Thus (jSnj holds for n = 2k + \. In the same way, we can show that it holds for n = 2k + 2. 

□ 

Proposition 6.3 There exists a non-zero constant a{E, h) ^ such that 

uq{x) = a{E, h) uq{x). 

The connection coefficients c^{E, h) in y28\) are analytic in E near Eq > and behave as 

c+{E,h) = a{E,h){l + o{l)), CQ{E,h) = -ia{E,h){l + o{l)) 
uniformly for E near Eq > as h ^ 0. 

Proof: The first part is a direct consequence of (|34|) and the construction of uo{x). 

The second part is an Airy type connection formula at least at the level of the principal 
term. Let us review briefly how to derive this by the Wronskian formulas (|21l) and ()22() . 
The coefficients c^, Cq are given by 

''wK+,Uo)' '° ''>V«,no-)' 

Since all the solutions involved just differ in the choice of the amplitude base point, we 
can apply formulas (|21|) and (|22|) to get 

>V«,Uo ) = 2i<™>^(:E;-;ro,:E+), W(no,no ) = 2i tq, 0). 

The Wronskian ^{uq^Uq) is more delicate, since there is a branch cut between the origin 
and x^, see Figure 2. Hence, Uq should be considered on the other Riemann surface. Let us 
denote by x the point on the Riemann surface continued from x to the same point turning 
around tq by the angle — 27r. Since gj^i^xY/"^ = —g^{xY^'^ and (7+(x)^/'* = igj^{xY^^, we 
have 

z(x,ro) = -z(x,ro), H{x) = -iH{x), T+{x) = iT_{x), 
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and for the series summands we get Wn,±{x) = w„^=p(x). Consequently, we have 

(x+ , ro , 0) = iu- (x+ ; ro , 0) , 

which yields 

W(no,^^^) = -2w;r°(x+;ro,0). 

On the other hand, we know that 

<-(xr;ro;x+) = l + 0(/i), 

because we can take a path from to xj' , along which x Rez(x) increases and which 
passes far away to the right from the turning point tq. Hence 

c+ = a<^'^(x-;ro; 0) (1 + 0(/i)) , Cq = ro; 0) (1 + 0(/i)) , 

and it is enough for the proof to show 

(42) lim wT"" {x7 ; ro ; 0) = lim wT"" (xt ; ro ; 0) = 1 . 

Let E>0. RecaU that w^^^f^ = E^=o(-^o ° )"(!)• Hence if we write x+ = iR, then 
ro;0) = = 1 + (/o o /i)«7)U.=,H 

where X{o,er/u) is the characteristic function of the interval (0, ER/u). The integrand func- 
tion is dominated by an integrable function: indeed, Lemma 16.21 gives for r G {0,ER/i') 

/ \ 2n / \ 2n 

n>0 n>0 ^ / 

|/i«™o'')(^'^^'/^)l < KcoshK, 
with = 1 + 3i?7^> and 

< X(0,i?ii/.)(O 1 + ^2(1 + ^2^2/^3)2 ^ 

On the other hand, the integrand function tends to as /i tends to (i.e. — > 0) for any 
fixed r > 0, since h{w'^'^f^)\x=o = 0. By Lebesgue's dominated convergence theorem, we 
obtain the second identity of (|42() for E > and by analyticity for E near £"0 > 0. The 
first identity of ()42p is obtained just in the same way. □ 

7. Connection formula near the critical point 

In this section, we compute the matrix T2, which connects uf{x) = u^{x]r\^xf) with 
uf{x) = u^{x; r2,xf), 

(n+(x),'u-(x)) = {u^{x),u~{x))T2{E,h), 

see (jnOJ. We will show 

Proposition 7.1 For all E near Eq > 0, the transfer matrix T2 is of the form 

t{E,h) s{E,h) \ 
-s{E,h) -t{E,h) J ' 



(43) 



T2{E,h) 
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and the asymptotics oft{E, h) and s{E, h) for /i — > are given by 

t{E,h) = D E-^'^ e-'^'^ + 0{h\\nh\), s{E,h) =i + 0{h). 

The proof of Proposition 17.11 is given in Section f7. 51 
Remark 7.2 From the symmetry properties 

(44) u+ = i ( ^ i ) '^r = -M ? n ) 



^ 1 7 ^' ' I V 1 . 

we know a priori that T2 is indeed of the claimed form (|43j) . Hence, it remains to prove 
the asymptotic behaviour of t{E, h) and s{E, h) as h ^ 0. O 



The exact WKB method is not enough to compute the /i-asymptotics of T2{E,h), 
because the two turning points ri and r2 tend to the same point x = \fE as /i — > 0. In 
other words: The Wronskian of two exact WKB solutions, say and u~, is given in 
terms of w^j"^" computed along a path from to passing between r\ and r2, but the 
/i-asymptotic formula for vf^^^ (Proposition 13. 3|) fails to hold because of the singularity 
of the function H at r\ and r2- 

Hence, one resorts to a microlocal study of the equation P^u = Eu near the point 
(x,^) = {^/E,^) for E > 0, where ^ is the dual variable of x. The equation P^u = Eu 
is reduced to a simple microlocal normal form Qw = (see Section [7.1|) . whose solutions 
are well studied. From these solutions we obtain two basis sets of microlocal solutions 
(/+, /~), {g^,g^) of Pi,u = Eu, which are related via a constant matrix R: 

(/+,/-)/? =(r,r) 

with 




see Section 17.21 The exact WKB solutions ti^ , uf are expressed in terms of these basis 
sets by 

{u+,uY) = {f+J-)Ai = {t,r)Bi, 

where the constant matrices Ai^r fmd Bi ^ satisfy 

(45) Ai = RBi, Ar = RBr. 
Then the matrix T2 is given by 

(46) T2 = A-^Ai = B-^Bi. 

Hence, the /i-asymptotics of t{E,h) and s{E,h) can be obtained from the study of the 
connection matrices A; ,, and Bi ^ (see Sections 17.31 and l7.4p . 
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7.1. Normal form. We now transform the equation {P,j — E)u = near (x,^) = {^/E,S,), 
> 0, to a simple microfocal normal form Qw = 0. 

Theorem 7.3 Let E > and u{x,h) be a solution of {Pi, — E)u = 0. Let V he the 
metaplectic operator associated with the j-rotation in phase space 

There exists a locally diffeomorphic change of coordinates x <P{x) = y, x > 0, with 
4>{y/E) = and a matrix-valued C°° -symbol M{y, h) = ld + 0{h) such that for any cut-off 
function x £ C^(M) identically equal to 1 in an interval around y = 

w{y, h) = V{x{y)M{y, h)u{(p'\y), h)) 

satisfies Qw = r, where 



Q 



7 = 'j{E,h) is a constant with -f{E,h) = -^E~^/'^h + 0{h^), and r{y,h) = 0{h°°) 
uniformly in an interval around y = together with all its derivatives. 

Remark 7.4 The terminology metaplectic operator will be recalled in Appendix El ^ 

Proof of Theorem\7^ We proceed in three steps to reduce the equation P,^u = Eu, that 
is 

x"^ — E V jx 



(47) hD^u{x) = A{x)u{x), A{x) 



-v/x E — X 



First step. One turns the quadratic diagonal entries of A{x) into linear ones. Let 
y = 4>{x) with 

1 /2 

c^{x) = {x - Ve) (^lix - Ve) + 2Ve'^ . 

In the complex right half-plane, the function </>(x) is a biholomorphic map with (j)[y/E) = 
and (/)(x)(/)'(x) = x'^ — E. The function 



(48) v(y) = HH^)) 



{1{x-./E) + 2VeY'^ 



x{x + y/E) 

is analytic in a neighborhood of y = 0, satisfying ■0(0) = and ip{(l){x))(j)' {x) = 1/x. 

Moreover, if u{x) satisfies (|17j) . then v{y) = v{4){x)) = u{x) satisfies 



Second step. Recall that v = vh with v G N — |. The second step makes the off-diagonal 
entries constant modulo 0{h^) by a change of the unknown function 

w{y,h) = M{y,h)v{y,h). 

Lemma ID . 1 1 constructs a matrix-valued C°°-symbol M(y, h) = ld + 0{h) such that w(y, h) 
satisfies 

(49) (^\^ ^^J^^w{y,h)=r{y,h)w{y,h) 
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where 7 = '^/'^h + 0{h?) and r(7/, h) = 0{h°°) uniformly in an interval around y = 

together with all its derivatives. 

Third step. Multiplying a cut off function x and then operating the metaplectic operator 
V from the left to equation (|49j) . we obtain by Lemma ID. 21 

Qw{y, h) = ix{y)r{y, h)w{y, h) - ihx'{y)w{y, h)) . 

The right hand side is of 0{h°°) uniformly in interval around y = together with its all 
derivatives. □ 

Remark 7.5 The proof of Theorem 17.31 shows, that 7, (p, and M depend analytically 
on E for E £ C near some Eo > 0. O 

Remark 7.6 In Appendix|E|we prove, that the M„(y) with M{y, h) ~ J2n=o Mn{y)h'^ 
have an analytic resummation, which is a Gevrey symbol of index 2. Hence, the second 
scale y/h, on which the two-scale Wigner measures in (TTj are based, reappears naturally 
in our normal form transformation. O 



7.2. Solutions of the normal form. Here we compute the solutions of the normal form. 
The equation 

Qw =0, w = *{WI,W2) 

is equivalent to 

(50) Wl = ■7=-W2 , -yW2 = -^W2 . 

V2y I 2 

This is a well-studied saddle point problem (see for example ^5]) [Z]i or section 5 in (29]). 
The system ()5U() has two maximal solutions 

with 

/f(y) = --^X(o,oo)(±y)|y|^i^i\ 
ifyv) = X(o,oo)(±y) lyl^'"^'' , 

where X(o,oo) is the characteristic function of the interval (0,oo). Moreover, we have two 
additional solutions 

(51) 5^(y)=(j \ ^CTd^{y\ 

where Thu(r\) = ^ e~^y^^^u{y)dy is the /i-Fourier transform and Cu = tZ is the 
complex conjugate. Indeed, we have the following identity: 

J I ) CJ^hQ = -Q ( J J ) C-^^- 
Since CJ^h = ^h^^: also have 

(52) f^{y)=(l l)cj-h9^{y). 
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These four solutions are linearly dependent. 

Proposition 7.7 The solutions and of Qw = are connected by 
(53) {g+,g') = {f+,r)R, R 

where 

P = --^r(l-4|7|^)exp(-|7|^), 



p q 

-q -p 



/l2 2hl'''l • r, TT r, 

— ^r(i-4h|>xp(--|,|=). 



Proof: The proof is just the one of Proposition 5.5 in 29 . We check g^ = pf^ — qf^ . 
One writes 

1 

gtiy) = CJ^hgiy) = hm -== / e^fe+-)V/^ T^-^kl^P^r? 

e-»0+ yj 1-Kh Jo 

and substitutes + ie)r]/h = —t to obtain 

where a(e) is the image of the interval (0, oo) under the map rj ^ {e—iy)rj/h. Using Euler's 
integral form of the gamma function T(z) = t^~^e~^dt, Re(z) > 0, and Cauchy's 
integral formula, one gets 

, l/2--L|7|2 

gtiy) = J- (-%)^l^l^-^r(l - 4|7|2) =p/+(y) _ qf-(y). 

□ 

Remark 7.8 In view of the relations H51|) and H52|) . the matrix R satisfies RR = Id, 
that is — = 1 and pq = pq. <> 

7.3. Frequency sets of the microlocal and WKB solutions. Let us study the fre- 
quency set of the microlocal solutions g^ and the exact WKB solutions u^^. 

First, the frequency set of the microlocal solutions are subsets of My x as follows: 

FS(/f) C {y = 0}U{r? = 0,±y>0}, 

FS{gt) C {?? = 0}U{y = 0,±?7>0}. 
Second, let a^^ C M^^ x be the Lagrangian manifolds defined by 
(54) at = {ro<x<ri,^ = ±{E - x^)}, at = {rs < x, ^ = ±(a:' - E)}. 



Since z'{x) = -y/i/^ — x'^{E — x'^y/x, one has by Lemma IC. 41 

FS('uf ) n {ro < X < n} C , FS('u^) n {x > C 



Now we transform to the normal form of Theorem 17.31 that is operate to the exact 
Nu{y,h) = V{x{y)M{y,h)u{r\y),h)) . 



WKB solutions uf^, 
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Since ^ are solutions of the equation {P^ — E)u = 0, we see that for any A; > 

D'yQiNnf^^) = 0{hn. 

Hence, Nuf^ are microlocal sohitions of Qw = near (0,0). Since the vector space of 
such solutions is two-dimensional, see Proposition 17 of [7], there exist matrix- valued 
C°°-symbols 

(55) A = (a;.fc), Ar = {a]^), A = (/3}fc), = [P]^), 

such that microlocally near (0, 0) 

{Nu+,NuJ) = {f^,r)Ai = {g^,g-)Bi, 
{Nu+,Nu;) = (/+,r)A, = {g+,g~)Br. 

Returning back to the (x,^) variables, i.e. operating A^^^ from the left, we have microlo- 
cally near (\/^, 0), 

(n+u-) = {f+J-)Ai = {g+rg-)Bu 
{ut,u-) = (/+,/- )A = {t,r)Br, 

where 

(56) f^ = N-\xf^), ~g^ = N~\x9^) 

and X £ C^(]R) a cut-off function, which is identically equal to 1 near y = and satisfies 
supp(x) C supp(x). 

By Lemma Id. 2| we have 

FS(/±) = k-1(a..FS(x/^)), FS{~g^) = (ac. FS(x<?^)) , 

where 

is the inverse of the canonical transformation K(^(x, ^) = {(j){x) , / (p' (x)) associated with (p. 
Clearly, 

K|FS(/f) C {rj = -y}U{rj = y,±y>0}, 

K^FS{gf) C {r, = y}U{i] = -y,Ty>0}- 

Since (/)(x)0'(x) = x'^ — E, this yields 

FS(/+) nU Cal^ U CJ+ U a~, FS(/") nU Ca+U a' U cr+, 

FS(5+) nU C U CJ+ U a~, FS{g') DU C cr^ U cjf U cr^ 

with U = {ro < X < ri or a; > r2}, see Figure 3. is neither contained in FS(nj^) nor in 
FS(/''"), while it is contained in FS(/~). Hence, = 0. Analogously one obtains 

(57) a{, = a^2 = /^L = = 0. 
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FSif-) FS(I+) 



FS (g- ) FS if-") 

Figure 3. Frequency sets of the microlocal solutions and defined in (|56() . 

7.4. Connection between microlocal and exact WKB solutions. We now compute 
the remaining coefficients of the matrices A; ^ and Bi r connecting the microlocal solutions 
with the WKB solutions u^, uf of P^u = Eu, E > 0. 

As a first step, the stationary phase method gives the following formulae for the mi- 
crolocal solutions, whose proof is to be found in Appendix IFl 

Lemma 7.9 For the microlocal solutions of PyU = Eu, E > 0, defined in l^5b]) . one 
has microlocally near the Lagrangian manifolds a'^ and defined in i)54l) 

(58) r{x) = e*"/« 2i/V"(^'"i)/'^ J ^ (1 + 0(/i)) 2ieara+, 

(59) /+(x) = e^"/8 2i/V^(^'"2)/fc ^ (1 + near a". 

It remains to connect e='=^(^'^i)/'^ and ^^^(^''^'2)/^ ^j^g exact WKB solutions uf^ and uf. 
Proposition 7.10 With the notation of Lemma 1 7. M we iiave microlocally 

= k^f^ near , = k^g~^ near , 

uj' = k^g^ near , = k^ near a~ , 

where 

k+ = -23/^6^8 (1 + 0{h)) , k; = _23/4e-i3-/8 (^^ ^ Q(^f^^^ ^ 

and 

— ji I I — li 
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Proof: To prove these relations, say the first one, it is enough to calculate the asymptotic 
behavior of and u'^ near a'^ . First, recall the exact WKB formula 

nf{x) = e±^(^'^^^)/'^r±(z(x;n)) (^|dd[^)^) = e^^^''^^^^ T^{z{x;n)) Q (l + 0{h)) 
with 

-L±[z) - ^ ^^(2) -H-^iz) ± iH{z) ) 

and 

Let X E]ro,ri[ be fixed. With the branches of the fourth root chosen in Section |5 one 
has H{z{x)) = e-^'^l^ + 0{K) as h ^ 0. Hence, 

and therefore 

(60) n+(x) = -2e^"/^e+^(^'^i)/'' J ^ (l + 0(/i)), 

n-(x) = +2e^'^/^e-"(^;"i)/'* J ^ (1 + 0(/i)) . 

Comparing (|En|l and we immediately obtain 

k+ = -2=^/V^/S(l + 0(/i)). 

Similarly, for fixed x in the interval ]r2,oo[, H(z{x)) = e'^*/^ + 0{h) as /i ^ 0, and 
utix) = +2e-^-/4e+^(-''^2)//. ( J ) (1 + 0{h)) , 

(61) n;(x) = -2e-^"/ V"^^''^^)/'' ° ^ (1 + 0(/i)) . 
Comparing (|HT|) and we immediately get 

fe- = _23/4e-i3-/8 (1 + . 

For the computation of and we use some symmetry properties. Recall that 
and 

= .^-^c, = vc, v-'T-^ = V, 

where the last identity comes from J^f^^ being the metaplectic operator of the transforma- 
tion (x,i^) I— > (^, — x), up to a normalizing constant. Hence, 
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and 
(62) 

On the other hand, 

u'^ = i 

which together with H62|) yields 



1 

1 



1 \ — 

1 T'- ' 



1 \^ 



1 



ki — iki , 



□ 



Proposition 17.101 yields for the connection matrices Ai^^ and Bi^r defined in H55|) 

0^21 = k'^ ) Oii2 = kj. , (322 — ^1 , (ill ~ • 

Combining these with the knowledge on the vanishing matrix elements ()57() and the rela- 
tions Ai = RBi, Ar = RBr, one obtains 

f k+ \ 

Br = 



(63) 



B, 



k7 



7.5. Computation of T2{E, h). Finally, we prove Proposition 17. ll calculating the asymp- 
totic behavior of t{E, h) and s{E, h) as , 

Let ^ > 0. By (gni) and 



(64) 



T2 = B^ ^Bi 



> 0. 






-Ph. \ 


q kr 


q kr 


p k^ 


q^ - p2 fc- 


q kr 


q kV / 



Remark 7.11 The identity ()64() is consistent with T2 



t s 
-s -t 



smce 



k; 



k; 



ky J kf Y k^' / kf 

which is checked by direct calculation. O 

By Proposition 17. 71 and Proposition 17. in( we then get 



p\ p 1 P^ — q^ 



Q q 



q kr 

Since 7 = -^E-^'^h + 0{h'^) and 



^/^r(l-4|7|2)eW4 



{l + 0{h)). 



sin(t|7P) ^^^^'^ 



l + 0(/i). 



we have 
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and 



' = -7:1^ = -' «^p(t ItI') (1 + o{h)) = -i + o{h). 



Since all the terms involved depend analytically on E for E £ C near Eq > 0, see 
Remark 17. 5| we have proven Proposition l7.il 

8. Proof of the main results 

In this section, we compute the Bohr-Sommerfeld type quantization condition of The- 
orem ^21 and derive the semiclassical distribution of resonances given in Theorem 11.31 

8.1. Quantization condition. Recall that Proposition 14. II gives the quantization condi- 
tion of resonances as c+(i?, h) = 0, where c^{E, h) is the product of three transfer matrices 
Ti , T2, T3 and the connection coefficients Cq , 

see identity ()32() . On the other hand, we have calculated the following asymptotics: 



^ (_^;^(,^ -+_^(^)), .^-;^e-W..^-3/^ + 0(.,ln.,), 



1 + 0(1) 

-i + o{l) 

see (jnSl, Proposition 17. II Proposition 15.31 and Proposition 16.31 respectively. Then, 

c+ = 2ae-^'/^e^^^/''{l + 0{h)) e'^«i/'^(l + o(l)) - e-^«i/''(l + o(l))) +0{e-^/^). 
Hence, c~^{E, h) = ii and only if 

(65) ^ V e-'-l^ E-^l^ e^^oi (^''^^/'^ + 1 = o(l) 

as /i — > 0, which proves Theorem 11.21 

8.2. Distribution of resonances. We now study the asymptotic behavior of the function 
Soi{E, h) as h ^ 0. Recall that 



1 _^2(^ 



r 



2\2 



Soi {E,h)= / ^ ^ dr. 



r 



If = P/i > is sufficiently small and E > 0, then i/^ — r'^{E — r^)^ > for r G (ro,ri), 
and the square root in the formula for Sqi {E, h) is taken to be positive. Substituting 
y = r'^ /E, one gets 



with uo = Tq/E and yi = rf/E. 
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Uo and yi are zeros of the cubic polynomial y{l — y)'^ — fi^ with fi = uE'~^/'^. If > 
is small and positive, then y(l — y)^ — ji^ has three zeros < yo(A*) < < 1 < y2(/^) 

with 2/0 (m) ~^ f-iid j/i,2(/^) — > 1 as /Lt ^ 0. We define 



where the square root is taken to be positive for < )U ^ 1. Since 



we study the asymptotic behavior of the function /(^u) as — > 0. For this, we have to 
understand the /U-dependance of ?/o(Ai) and yi{l^)- When /j? turns around once in the 
positive sense (i.e. /j, becomes e'^*/x), then yo{f^) turns around in the positive sense and 
yi(/i) and y2{p) exchange their position turning half around 1 in the positive sense. As a 
consequence, taking the branch into account, 

(66) /(e» = /(/x) + it;(M) + r(/x) 

with 

. f Vm2_^(i_^)2 ry.i,) ^^2_^(i_^)2 

^(y") = -M ^ c^y^ ^(m) = / ^ dy, 

where Tq is a contour around 0. These functions have the following properties: 
Lemma 8.1 T'(/x) is a holomorphic function of fi^ at fi = 0, and 

i?(/x) = -TT//, r(M) = ^(i + o(/.2)). 

In particular, i?(e''V) = -R{l^) and r(e''V) = ^^C/^)- 



Proof: By the residue theorem, R{p) = t^\/ fJ-"^ = —nfi. For the study of T(/i), wc move 
by the locally biholomorphic change of variables v = ^Jy{y — 1) from a neighborhood of 
y = 1 to a neighborhood of u = 0, 



T{^l) = -i r ^/fJ,■^-v^ f{v)dv 



where f{v) = {2y{v)) ^-^yiv) is holomorphic in a neighborhood of f = and satisifies 
f{v) = ^ + f'{0)v + 0(i;2). Hence, 

T{fx)=if/J ^/l - fifxw) dw = if? (^^ J VT^^ dw + Oifi^)^ , 

since J^^ \/l — w'^ w dw = 0. □ 

Proposition 8.2 /(/it) is ramified at /^ = and satisfies 

/(M) = | + fM + 0(/x2|ln;,|) (/x^O). 
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Proof: By Lemma l8. 11 we have from (|66|) 

/(e2-V) = /(e"V) + ^(e"V) + r(e"V) = a(/")+^(^)+r(^))-i?(^)+r(^) 
= /(/i) + 2r(/i). 

Since ln(e^'^*//) = In fi + 27ri, this means that the function 

^(^) =/(/.)- 4 r(^) In ^ 

is single- valued around fi = 0. Moreover, since T{fi) behaves quadratically in fi near = 0, 
A{iJ,) is holomorphic near fj, = with 

A{fx) ^ A(0) = 7(0) = / (1 - = |. 

Differentiating equation (|(i6j) . one gets A'(0) = /'(O) = ^. Hence, 

= A(0)+yl'(0)^+i-T(/x)ln/x + O(/x2) = | + i.r(^)ln/i + f// + O(//2) 
= | + f^ + 0(/i2|ln/i|). 

□ 

Proof of Theorem \l.,H The quantization condition (|65|) is satisfied, if and only if there 
exists an integer A; G Z such that 

(67) 2Soi{E,h) -^hlnE^/^ + ^hln{^9^h) = (2A: + |) m/i + o(/i). 
Setting A = E^^"^, Proposition 18 . 21 implies 

Soi{E,h) = iXI{9h/X) = liX + ^i9h + 0{h^\lnh\), 

and the quantization condition 1)671) becomes 

|a-i/ilnA + i/iln(fp2/i) = (2A;-z; + |)i7r/i + o(/i). 

Writing A = Ai + 1X2 with Ai, A2 G M, the real and imaginary part of the above condition 
read as 

(68) |Ai-i/iargA = {2k - i) + I) irh + o{h) , 

(69) -|A2 - i/iln|A| + i/iln(fi>2/i) = o(/i). 

Now, we assume that a < Ai < 6 and A2 = o(l) as h ^ 0. Then, 

argA = arctan(A2/Ai) = o(l), ln|A| = InAi + | ln(l + A2/A?) = InAi + o(l). 
Setting Xku = ff (8A: — 4P + 5), euqations (|H5|) and become 

Ai = %{8k-4i) + 5)h + o{h) = Xki>h + o{h), 

A2 = -i(/^lni-Mn^)+o(/.). 

□ 
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Appendix A. Energy surfaces 

Let £^ G C and u e h{Z + l). The zeros of the function C+ ^ C, r ^ {E - r'^ f - v"^ / r"^ 
are the roots of the sixth order polynomial — 2Er'^ + E'^r'^ — v"^ in r, which lie in the 
right half-plane C"*" = {r G C; Re(r) > 0}. This polynomial has at most three different 
roots ro,ri,r2 S C"*" in the right half-plane, whose squares are the roots xq,xi^X2 G C of 
the cubic polynomial — 2Ex'^ + E'^x — in x. \i E ^ M, then an easy criterion for 
real-valuedness of the roots xq^x\^X2 is the sign check of the polynomial discriminant 

= (XO - Xif (Xo - X2f (Xi - X2f = (iE^ - 271/2) . 

The three roots are real if and only if D3 > 0, that is iff < 4E'^/27. The roots are real 
and distinct, if and only if u"^ < 4i?^/27. 
From Cardano's formula 

xq = 2E/3 + S+ + S^, 

xi,2 = 2E/3-{S+ + S-)/2±iV^{S+-S^)/2 

with S± = ^-E^/27 + 1/2/4 ± ^^Ds/lOS one learns xq ^ and xi,2 ^ E as h ^ 0. 

If one denotes by G the square root of G C, which lies in the right half-plane, 
then ro — > and ri^2 — > VE as /i — > 0, see also Figure 4. 



Energy surfaces for E=2 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 



Figure 4. The energy surfaces {(r, ;o) G M+ x M; = [E - r^Y - v"^ /r"^} 
for E = 2 and different values of v G {^ij?!}- 

Appendix B. Spectrum of P+ = -/i^A -|- |x| 

The Schrodinger operator P"*" has a locally bounded positive potential, which increases 
to infinity as |x| 00. Hence, is essentially self-adjoint on C^(M2) and has purely 
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discrete spectrum. We are looking for eigenvalues E £]a, b[ in a bounded positive interval 
]a, b[c M"*". In polar coordinates x = r(cos 9, sinO), r > 0, G T the differential expression 
—h?IS.x + |x| reads as 

Hence, E \s a. formal solution of the eigenvalue problem (P"*" — E)'il: = if and only if there 
exists / G N U {0} such that 

, 2 / d^ Id f\ , , , , , 

\-^2+-^,- -2) Mr) + {r- E)wi{r) = 0. 

Substituting wi{r) = r~^/^n^(r), this is equivalent to 

(70) -h" + ui{r) + (r - E)ui{r) = 0. 

The ordinary differential equation (|70|) has r = as a regular singular point with 
exponents ^ ± ^, while r = oo is an irregular singular point of rank two. There are two 
linearly independent solutions, which behave as r^^"^ and r^/^ Inr (for / = 0) or r^/'^+^ and 
7.1/2-^ (Jq^ I near r = 0, see for example Chapter 5.5 in |28j . 

Looking at the irregular singularity at infinity, there is a fundamental system n^(r), 
which behaves as {r—E)~^/^ exp(ib(r — i?)^/^//i) near r = 00. Hence, a necessary condition 
for E being an eigenvalue of reads as follows. 

Lemma B.l E G o-disc(-P^) if and only if there exists I G NU{0} such that >V(n°, u^) = 
0, where W{u^,u^) = {u^)' — (u^)' is the Wronskian of the two solutions to ifTDj) . 

n-(r) ~ (r-i?)-i/4e-(r-i?)^/V^^ r^oc. 



Let us have a closer look at the potential of equation (|7U|) 

Q(r) = r-E + h^{l^ - l/4:)/r^. 

The origin r = is a double pole and the coefficient h'^{l'^ — 1/4) is positive for / > 
and is negative for / = 0. This difference is crucial from the WKB point of view, since 
the level curves of Re(J y^Qir) dr) in the complex r-plane are closed curves enclosing the 
origin for / > and radial curves emanating the origin for 1 = 0. The WKB construction 
at a double pole for / > has been treated in [14^, and we here restrict ourselves to the 
application of this results. 

We assume / > 0. For h > sufficiently small, the potential Q has three distinct simple 
turning points: they are zeros of the cubic polynomial — Er^ + h'^{P — 1/4). Two of 
them are at a distance 0(h) from r = 0, while the third is at a distance 0(h) from r = E. 
Since the polynomial discriminant 

h\l'-l)[^E-^-i^(l'-\))>0 

is positive, all three turning points are real. The first one is negative, while the other 
two are positive. We denote these two by 01,02 > 0. The strategy for characterizing the 
quantization condition yV(u^,u^) = is the following. We connect 
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(1) uP with two exact WKB solutions built with the Langer-modified potential 

Q{r)=r-E + h^l^/r^ 



(3) with an exact WKB solution u^. defined for r > 02, which is collinear to u^. 



defined for r G]ai,a2[ with the phase base point at ai, 
(2) with two exact WKB solutions uf^ defined for r G]ai,a2[ with the phase base 
point at 02, 

We start with the connection near the origin. Proposition 11 in ^Ij proves existence of 
a non-zero constant a = a{E,h) 7^ such that 

u° = a (1 + 0(1)) 5+ + a (1 + 0(1)) u- 

as /i ^ 0, where 

5±(r) = g(r)-i/4exp(^±-^ Qis^/'^ds^ W^{r;r^) . 

5i > denotes the first positive turning point of Q, while rf^ £ C are two suitably chosen 
/i-independent points in the complex plane such that Re{r^) €^]ai,a2[, ±lm{r^) > 0, and 

Next, we connect near the turning point. The WKB solutions u^{r) = u^{r;a2,r^) 
are of the form 

with W^{r^;r^) = 1. Connecting and uf, 

we evaluate the Wronskians, which are constant with respect to r, near the points r = r^. 
By this choice we can exploit the fact that Q{r) = Q{r) + 0{h?) as h ^ uniformly 
in neighborhoods of r = r^. Also the turning points differ by a term of order /i^, 5i = 
ai + 0(/i^), and one derives 

5±(r) = u^{r;ai,r^) (1 + 0(/i)) 

as /i ^ uniformly in neighborhoods of r = r^. Hence for the connection formula, we 
just have to account for solutions with different phase base points ai, 02- By Proposition 
2.5 in Uni, one has 



W{u+,ui) = -l{l + 0{h)), 

W{u^,uf) = e±^i2/'^(l + 0(/i)) 
where 5 > does not depend on h, and 



Si2 = Si2iE,h) = / Q(r)i/2dr 

J ai 



is the action integral between ai and 02- We note, that the square root is taken such that 
5i2 is purely imaginary with positive imaginary part. Collecting the previous formulae, 
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one obtains 



/e^^^/'^ (1 + 0{h)) 0(6-5^) 
[u ,u )-[^Ui,u,)y o(e-<5/^) (i + 0{h)) 

and 

(71) = ae^'^l^{l + o{l))u+ + ae-^^^l^{i + o{l))uJ (/i ^ 0). 

Finally, we connect to infinity. The solution u^{r) = (r; a2,rr) is defined as 



u,. (r) 



with G M such that > 02 and W~ {rr;rr) = 1. Since 

Q{r) ^ r — E, r — > 00 
u~ is cohinear to u'^, and W(ti'^, u^) = if and only if y\^{vP ,u~) = 0. By equation ((7T|) 



it is enough to compute the Wronskians W{uf,Uj. ). Choosing a branch cut for (5(r)^/^ 



near 02 as in Figure 2 of one gets 

>V(n+ n,-) = -f (1 + 0(/i)) , W{uJ,u-) = | (1 + 0(/i)) , 

since {r) = — i?i"'"(r; 02, S;~), where is the point on the other Riemann surface. 
Hence, we have proven 

Proposition B.2 Assume I > 0. Then, 

>V(n°,u-) = =^ e^^i^^^''^)/'^ + 1 = 0(1), h^O, 

with 

(72) Si2{E,h) = n ^r-E + h?{P -\)/r^ dr. 

What remains, is a study of the asymptotic behavior of the action integral 6*12 (-E, h) as 
/i ^ 0. Substituting r = Ey in (|72j) . we rewrite 

5i2 (i?, /l) = 1 (^) , /X = ^^^3~,^^^ /I, 

with 



'2/1 (m) y 

where yiip) and 1/2 (a*) are the first and second positive zero of the cubic polynomial 
y'^{l—y)—/j?. For /i > small enough, the polynomial has three zeros yo(/^) < < yi(/u) < 
y2{fJ') < 1 with yo{fJ-),yi{lJ') and y2(Ai) — > 1 as /x — > 0. When turns half around zero 
in the positive sense, that is when /x turns to e"^/i, then yo{fJ-) and yi(/x) exchange their 
position, while y2{fJ') turns around 1. Hence, /"'"(e"^//) = + R^ifJ-) + T~^{fi) with 

v5EZIIEi),„ 



and T~^{^i) = /p^ y^y2(l — y) — fi'^/ydy = 0, where Fi is a contour around 1. Substituting 
V = y-^/l — y, we rewrite 



R+{f,)=-i ^^j? -v^ ^ dv, 

J fir\—fi ^ 
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where f{v) = l{v^ — \y(v)^') is holomorphic in a neighborhood of f = and satisfies 
/(v) = 1 + i„ + /"(0)w2 + 0(v^). Then, 

Jlrx-l ^ 

as /i ^ 0, since j-^^_^^/T^^iv^ dw equals — ivr, — 1,0 for j = —1,0,1, respectively. 
Hence, 

/+(e-/i) = - vr^ + f ^ 0(^2)) ^ 

and Proposition 18.21 yields 

/+(M) = | + f^ + 0(/i'|lnH) (/U-0). 

Remark B.3 With these asymptotics, the quantization condition of Proposition IB. 21 
can be rephrased as follows. There exist an integer A; G Z and a natural number / G N 
such that Si2{E, h) = nr{k + + o{h) as /i ^ or 

E^/^ = Itt (2k + I - y^P^Tjtj h + oih) (h^O). 

O 

The study of the resonant set of the full operator P and the spectrum of the upper level 
operator P'^ have revealed a parallel structure: A resonance £" of P is characterized by a 
Bohr-Sommerfeld type condition 

where t is related to non-adiabatic transitions, and SQi{E,h) is an action integral. An 
eigenvalue A of P^ is characterized by a Bohr-Sommerfeld condition 

where 5*12 (A, /i) is an action integral. The integrals can be expressed as 
Soi{E,h) =iE^/^I{n), fi = hi)E-'^/^, PGN + i, 

where /(/i) and /^(/i) share the same /U-asymptotics for — > 0. 

Appendix C. Symbols, frequency sets, microlocal solutions 

Let us briefly recall the definition of the symbol classes used, as well as the notion of 
frequency set and microlocal solution. 

Definition C.l A smooth function a{x, h) on an open set J7 C M is called a C°°- 
symbol of order m G Z, if there exist smooth functions an{x) G C°°(f2,C) such that for all 
non- negative integers a G N and G N, 

N 

sup I d^ia{x, h)-Y, a„(a;)/i"+"^) | = 0(/i'^+^+^), h ^ 0. 

x&n 

When this holds, we write a{x,h) ~ Xl^o '^"(^)^"~'''"- 
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Definition C.2 An analytic function a(x, h) on an open set C C is called a Gevrey 
symbol of index 2, if there exist functions an{x) analytic in 0, such that for all compact 
subsets K C 0, and all large enough C > there exists 6 > with 

1 

sup I a{x, 



Definition C.3 Let u £ S' he a possibly /i-dependent distribution and (xc^o) G T*W. 
Then (xo,^o) ^ FS(u) if there exists a Schwartz function xo on phase space T*M with 
Xoi^o, ^o) = 1 such that for any N £ N 

Xo{x.hD^)u = 0{h^), h^O. 

FS{u) is called the frequency set of u, and u a microlocal solution of u = near (xo,Co)- 

In Sections 17.31 and [7^ we have used the following fact about the frequency set of WKB 
solutions: 

Lemma C.4 If u{x, h) = a{x^ h) eiqi{i(j){x) / h) where (f) G C^{Q.^ M) is a real phase and 
a is a C'^-symbol, then FS(u) C {{x,S,); ^ = d.j,(l){x)}. 



Appendix D. Proof of Theorem 17.31 

Lemma D.l Let P G N — ^, y ^ ip{y) the function defined in i)48|) with ^/'(O) = , 
E > 0, and v{y) = v{y, h) a solution of 

There exists a matrix-valued C°^-symbol M{y,h) = Id + 0{h), such that w(y,h) = 
M{y,h)v{y,h) satisfies 

(74) (^^^_" y ^^J^ ^ w{y, h) = r{y, h)w{y, h), 

where 7 = -^E~^/^h + 0{h'^) and r{y, h) = 0{h°°) uniformly in an interval around y = 
together with all its derivatives. 

Proof: We rewrite equation (f75|) as hDyv{y,h) = B{y,h)v{y,h) with 
and equation (|7i|) as {hDy — G{y, h)) w{y, h) = r{y, h)w{y, h) with 

In 



G{y,h) ~ Go(y) = Bo{y), ^ ( _° ) (n > 1). 



ra=0 

We are looking for 



00 



M(y, h)^Yl Mn{y)h\ Mn{y) = r'"]'^] (n > 0) 
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such that 



or equivalently for all n > 



hDyM = GM - MB 



(75) DyMn-l = (GjMn-j - MjBn-j) 

j=0 

with convention M_i = 0, i?„ = for n> 2. 

Equation ((7^ is satisfied for n = 0, if we take Mq^u) = Id, i. e. mo{y) = 1 and qo{y) = 0. 
Then, we have for n > 1 

Gn = -iM^_^-Y,GjMn-j + Mn-lBi + MnBo 
j=0 

/ n \ 

= -iM;_i - ^ G,Mn-j + Mn-iBi - 2y i '^^ j . 

Since V'(y) ^ for y G M, the previous equation is equivalent to 

n-l 

(76) 7n = -iq'n^iiy) -^jjuin-jiy) + i>'il^{y)mn-i{y) -2yqn{y), 

ra-l 

(77) = -im'^-i{y) -YuQu-jiy) - i>iJiy)qn~i{y)- 

i=i 

Let us start with n = 1. Substituting y = in (jTB)) yields 

71 = t'V'(O), 

and we automatically obtain 

Qiiy) = -77- (71 - J^i^iy)) , 

which is smooth (even analytic) near y = 0. Setting mi(0) = 0, equation (|77() gives 

rv 

fni{y) = i I (71^1 (y') + j>'ip{y')qi{y')) dy'. 

Jo 

If ^j,mj,qj with mj(0) = are determined for 1 < j < n, then (fTH]) yields 

7„+i = -ig^(O) 

and 

1 " 
qn+i{y) = -^(7n + «g[i(y) + - t'V'(y)"i„.(y)), 

^ i=i 

which is smooth (even analytic) near y = 0. Setting m„4.i(0) = 0, equation ((77j) then 
gives 

»j, n+l 

"i„+i(y) = i / (V 7j9n+2-j(y') + ^'^{y')(ln+i{y'))dy' . 
Jo 

Since mn{y) and qn{y), > 0, are smooth functions near y = 0, there exists a matrix- 
valued C°°-symbol M(y, /i) with M(y, /i) ~ Y^^^q M„(y)/i" such that CZH) holds. □ 
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Lemma D.2 Let K^{y,ri) = -^{y — V^V ~^ v) be the j-rotation in phase space T*R. 
The metaplectic operator V of the transpose hnear canonical transformation 

4 4 

satisfies 

V{hDy -y) = -V2yV , V{hDy + y) = V2 HDyV. 
Moreover, FS(Fn) = k_|FS(u) for u G S' . 

Proof: One applies Theorem 2.15 in jl2j for the hnear canonical transformation kil. □ 
Remark D.3 A formula for V in terms of oscillatory integrals is given by 

Jm. 

see Theorem 4.53 in |12| or Proposition 5.3 in [2^1 . O 

Appendix E. Gevrey symbols 
Lemma E.l The analytic functions mo{x) = 1, qo{x) = 0, 



qn{x) 



rrinix) 



1 / /"^ 1 



-mn-i{x)i'ip{x) + mn-i{x) D^{0) 



n > 1, 



i> Qn{t)i^{t)dt + qn+i^j{t)dt, n>l 



with 7i = i>'i/'(0) ^.nd 7„ = —iq'j^_i{0) for n > 2, which have been introduced in the proof 
of Lemma lD.it have analytic resummations q{x, h) and m{x, h) in some neighborhood of 
X = 0, which are Gevrey symbols of index 2. 



Proof: Let il. = {x £ C; \x\ < r} and = {x G C; |2;| <r — t} for < t < r. It is enough 
to show, that there exist positive constants D^, Dq, C > such that for all n G N 



(78) 



sup \qn{x)\ < DqC 



.(2n) 



2n 



t 



2n 



sup \mn{x)\ < DmC 



,(2n) 



2n 



f 



2n 



We use the following bounds: if f{x) is a holomorphic function in which satisfies for 
some positive constant M > 

M 



sup |/(x)| < -TT, 



then 



and 



sup 



f{u)du 



sup sup 



< 



< 



M 



sup 



f{u)du 



< 



M 



, ,,,, (k + 2)''+^M 
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We prove ((75|l by induction: 



n-l 



sup \qn{x)\ < sup \q'^_^{x)\ + ^ kj_i(0)| sup 



+2|P| sup 



sup |^(x)| 



n-l(2n) 



2n-l 



n-i(2ri-ir 



— sup|V'(x 

< DJl + D,r. + 2\D\r^ sup mx)\) C^~' ^^""^ 
V xen / 



f2r, . 

\2n 



where we used the inequality 

k\N - kf-"" < {^f {f)<k<N). 

Taking C large enough so that C > 1 + Dm + 2r|P|-^ sup^g^ IV'l^^)!; we obtain the first 
bound in ()78() . Similarly, one has 



(2n 



,2n 



sup |m„(x)| < sup|V'(a;)|i:'gC" 1^,2n-l 

71 



+ 2^i^,0 (ij ij iJgO (2n + l-2j)t2n 



n(2n) 



2n 



< (|i>|r sup|V(x)|+2D,)i?,C"l^. 



Taking Dm so that D„i > (|z>|rsup^gQ + '2Dq)Dg, we get the second bound in (f75)) . 

□ 

Appendix F. Proof of Lemma EH 

In this appendix, we prove one of the two formulae of Lemma 17.91 i.e. compute the 
asymptotic exapnsion of microlocally near . 

Recall that near y = 

f+{cp-\y)) = M{y,h)-W-^f+{y) 
where the inverse of the metaplectic operator is given by 

y-l f+{y) = e-/8(V2^/i)-l/2 f ^i{y^+z^^2^yz)l{2h) f+^,)dz, 



and that 



V2y 



X(o,oo)(y) lyhh'"^'" , f2 = X(o,oo)(y) \y\ 



Let us compute the asymptotic expansion of 

l-OO -1 

V~^f+{y) = e^"/8(^/27r/i)-i/2 / ^My,z)/h^^^ ^^y^^^ ^ _^y2 ^ ^2 _ ^^^^ ^ \^f\xiz). 

Jo 2 
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Since 7 = 0{h), the phase function ip{y,z) has two real critical points 



zt{y) = {y±Vy'-h\^)/V2, 

which are the roots of ^{y,z) = z — \/2y + The phase of the asymptotic expansion 
will be given by ip{y, zf{y)). Since we are on a~ , we can assume x > VE, i.e. y > and 
independently from h. Then, 

|7P 



and 



zt{y) = V2y + 0{h% .-(y) = J-^ + 0(/i2) 



ip{y,zt{y)) = Tly' + 0{h^) 



In view of Lemma IC. 41 this means that the critical points z^{y) and z^ (y) will contribute 
on a~ and a^, respectively. Hence, we have only to compute the contribution from z^{y). 

Since z+(y)) = 1 — ^2 = 1 + 0{h'^), z^{y) is a non-degenerate critical point, 

\y} 

and the stationary phase theorem (e.g. Proposition 5.2 in j^I) says that 

microlocally near a~ . Comparing iLp(^(^(x) ^ z^ {<\){x)y) = — + 0(/i^) with z{x;r2) = 
J^^ — t'^{E — t'^Y /tdt, we observe for x > r2 

-^z{x-r2) = -^1/2 _ x2^E - x2)2/x = -i{x^ - E) + 0{h^) 
ax 

since u"^ = 0{h?) and 0(x)0'(x) = x'^ — E. Because of z{r2;r2) = 0, one obtains 

-z{x; ra) = -^(^ixf - ^(p{r2f + 0{h^). 
If (/>(r2)^ = 0{h?), then i(/9((/>(x), = -z{x;r2) +0{h'^). Consequently, 

V-^f+{(l){x)) = e*"/s 2^/4 e-^(^''-2)/'^ (1 + 0{h)) 

and 

microlocally near o"", since M{(j){x), h) = Id + 0{h). 

It remains to prove 4>{r2)'^ = O^h?). The turning point r2 is a zero of the function 
X I— > jx^ — {E — x'^Y'. This is equivalent to det(a(r2,0)) = with 



the symbol of the operator A{x) — hD^- The first step of the normal form transformation 
of Theorem 17.31 reads on the symbol level as 

a(x,0 = (li'{x)a{(t){x),i/(t)'{x)) 

with 



-i'ijj{x) — X J 



44 



S. FUJIIE, C. LASSER, AND L. NEDELEC 



Since (t>'{r2) 7^ 0, one has det(a(r2,0)) = if and only if det(a((/)(r2), 0)) = 0. The rest of 
the normal form transformation is 

with 




the symbol of the normal form operator Q, and jj/i the Moyal product of semiclassical Weyl 
calculus. The /i-asymptotics of fth (e.g. Proposition 7.7 in 9 ) together with the linear 
^-dependance of a(x,^) yield 

MUaU M-^ = MaM-^ - ihW = MaM'^ + ©(/i^). 

Hence, a = V2M-'^{q o K|)M + ©(/i^) and 

det(a(r2,0)) = ^ det((7 o (ra, 0)) = ©(/i^) <^ ^(ra)^ = 0(/i^). 
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